function R = AIA_step1_R(I,Delta,r,c)
    F = length(I);
    R11 = 0;R12 = 0;R13 = 0;
    for f = 1:F
        delta = Delta(f);
        i = I{f,1}(r,c);
        R11 = i + R11;
        R12 = i*cos(delta) + R12;
        R13 = i*sin(delta) + R13;
    end
    R = [R11, R12, R13];
    R = R';
end